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Abstract. We present a simple framework to compute hyperbolic Voronoi 
diagrams of finite point sets as afhne diagrams. We prove that bisectors 
in Klein's non-conformal disk model are hyperplanes that can be inter- 
preted as power bisectors of Euclidean balls. Therefore our method sim- 
ply consists in computing an equivalent clipped power diagram followed 
by a mapping transformation depending on the selected representation 
of the hyperbolic space (e.g., Poincare conformal disk or upper-plane 
representations). We discuss on extensions of this approach to weighted 
and fc-order diagrams, and describe their dual triangulations. Finally, we 
consider two useful primitives on the hyperbolic Voronoi diagrams for de- 
signing tailored user interfaces of an image catalog browsing application 
in the hyperbolic disk: (1) finding nearest neighbors, and (2) computing 
smallest enclosing balls. 



1 Introduction and prior work 

The birth and spread of science was originally initiated by geometry, the science 
of (Earth) measurements. Euclid (300 BC) laid the foundations of geometry in 
his masterpiece book series Elements that still remain the basis of mathematics 
nowadays, two millcnia later. Euclidean geometry was assumed to be the only 
consistent geometry until the nineteenth century. Eventually, failing to prove 
Euclid's fifth postulate (known as the parallel postulate) from the others opened 
the door to abstract geometries. The first two abstract geometries that were 
historically called imaginary geometries are the hyperbolic geometry and the 
spherical geometry. The essential difference between Euclidean and these non- 
Euclidean geometries is the nature of parallel lines. While in Euclidean geometry 
there is a unique line passing through a given point and parallel to another 
line, there can be infinitely many in the hyperbolic geometry [1] and none in 
spherical/elliptical geometries. In this paper, we focus on characterizing and 
building efficiently the Voronoi diagram in hyperbolic geometry. 

The Voronoi diagram [2] of a finite point set partitions the space in cells 
according to proximity relations induced by a distance function d. The structures 
of Voronoi diagrams appear in various fields such as crystallography, astronomy 
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and biology. Let V = {pi, ...,p„} be a set of n <i-dimensional vector points of 
space X . Then the Voronoi diagram of V is defined by the collection of proximal 
regions vor(p i )'s, called Voronoi cells, such that vor(pj) = {x e X \ d(pi,x) < 
rf(p 3 ,x), Vp 3 e?}. 

The Voronoi diagram in hyperbolic geometry has already been partially stud- 
ied; Onishi and Takayama [3] investigated the hyperbolic Voronoi diagram con- 
struction in the Poincare upper half-plane. Onishi and Itoh [4] further extended 
the Voronoi diagram in Hadamard manifold that are simply connected complete 
manifold with non-positive curvature. Boissonnat et al. [5] considered the hyper- 
bolic Voronoi diagram in the Poincare conformal upper plane and the Poincare 
conformal disk. 3 Boissonnat and Yvinec described in their computational ge- 
ometry textbook the hyperbolic Voronoi diagram in the Poincare <i-dimcnsional 
half-space [2], pp. 449-454. They do not introduce explicitly the hyperbolic 
distance but rather show how to answer proximity queries (e.g., whether B or C 
is closer to A) using the notion of pencil of spheres. They deduced that the com- 
plexity of the hyperbolic Voronoi diagram is 0(nlogn + n^ Nilforoushan and 
Mohades [6] concentrated on the hyperbolic Voronoi diagram in the Poincare 
2D disk for which they report an incremental quadratic algorithm, and char- 
acterize geometrically the orthogonality of bisectors with geodesies. Note that 
by using the graphics processor unit (GPU), it is easy to rasterize interactively 
Voronoi diagrams [7] for any arbitrary distance function, including the distance 
of hyperbolic geometry. 

In this paper, we revisit these work by considering the Klein projective disk 
model. We showed in section 2 that bisectors are hypcrplancs, implying that the 
Voronoi diagram is affine and therefore can be easily constructed from an equiv- 
alent power diagram. We report the one-to-one conversion formula to change 
the representation of the hyperbolic geometry: Klein or Poincare disk models 
and Poincare upper half plane model. We further characterize the weighted and 
fc-order diagrams in the Klein disk model and explain the dual hyperbolic De- 
launay triangulation. In section 3, we present an image browsing application in 
the hyperbolic disk that requires two basic user interface selection operations 
that are solved by means of nearest neighbor queries and by finding the smallest 
enclosing ball in hyperbolic geometry. 

2 Hyperbolic Voronoi diagrams from power diagrams 
2.1 Klein's non-conformal disk model 

The Klein model (or Beltrami-Klein model [1]) uses the interior of a unit circle 
for fully representing the hyperbolic plane. The Klein model is also known as the 
projective disk model where lines are depicted by chords of the circle (e.g., line 
segments joining any two points of the circle). Although the Klein model offers 
a simple visualization of geodesies as line segments, it has the disadvantage of 
being a non-conformal representation: That is, angles of the hyperbolic plane are 

3 See Figure 6 and Figure 8 of [5] for illustrations. 
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not preserved, and get distorted in this model. (That explains why the Poincare 
conformal disk representation is often preferred.) The distance between any two 
points p and q of the Klein disk with Euclidean coordinates (||p|| < 1 and 
||q|| < 1) is computed as 

^k(p, q) = arccosh 1 (P'°t) n\ 

va-iiPipxi-iiqii 2 )' 

where arccosh a; = log (a; + \J x 2 — 1). 



2.2 Bisectors in Klein projective model are hyperplanes 

Bisectors are core primitives to characterize and compute Voronoi diagrams. 
The bisector B of two points p and q with respect to a given distance function 
D is defined as the loci of the points c that are at the same distance of to p 
and q: B : {c | D(p,c) = -D(q, c)}. Let us derive from the distance function of 
Eq. 1 the equation of the bisector in Klein's disc model of hyperbolic geometry: 
h K (p,c) = /i K (q, c). That is 

1 ~ (P, c) = 1 - (c, q) 

V(l-||p|| 2 )(l-||c|| 2 ) V(l-Hc|| 2 )(l-||q|| 2 )' 
(c, y/1 - ||p||2q- y/1 - ||q||2p) + y/l - ||q||2 ||p||2 = (2) 

This is a linear equation (a, c) + b = in c, with a = yfl — ||p|| 2 q— \J\ — ||q|| 2 P 
and b = yjl — ||q|| 2 — \J\ — ||p|| 2 - Thus the bisector is a hyperplane, and the 
Voronoi diagram is said affine [8]. 



2.3 Hyperbolic Voronoi diagrams as equivalent power diagrams 

The power distance of a point x to a Euclidean ball B — B(p,r) is defined as 
||p — x|| 2 — r 2 . Given n balls Bi — B(pi, n), i — 1, . . . , n, the power diagram (also 
called Laguerre diagram) of the Bi 's is defined as the minimization diagram of 
the corresponding n functions D^x) = ||pj — x|| 2 — r 2 . The power bisector of 
any two balls B{p il r i ) and B(pj,rj) is the radical hyperplane [8] of equation 

2(x,p, - Pi ) + Hp.II 2 - ||p, 1 1 2 + r) - r\ = 0. (3) 

Thus power diagrams are also affine diagrams. Interestingly, Aurcnham- 
mer [8] proved that all affine diagrams can be constructed from equivalent power 
diagrams. That is power diagrams are universal representations of affine dia- 
grams. Therefore any algorithm to build power diagrams can also build all affine 
Voronoi diagrams once the equivalence relations are explicitly given. Thus we 
need to report the mapping that associates to points in the Klein disk corre- 
sponding Euclidean balls so that their bisectors match. By identifying the bisec- 
tor equation Eq. 2 with that of Eq. 3, we find that point Pi is mapped to a ball 

of center Ci = — . Pi and radius r, = ,,} P u ^ H2 -> , 1 . Note that some 

2V1-llp.ll 2 4(i-|| Pi |p) ^A-IIp*!! 2 
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of the radii may be negative (imaginary numbers, e.g. ||p|| 2 = — > r = — 1 = i 2 ), 
i.e. represent imaginary balls. A careful analysis shows that the radius is negative 




— 1) ~ 0.9718... and positive otherwise. 



Theorem 1. The hyperbolic Voronoi diagram of n d- dimensional points can be 
obtained in the Klein disk model as an equivalent power diagram. The diagram 
has combinatorial complexity 0(n'^) and can be built in 0{n\ogn+n^^) time. 

Since power diagrams are defined on the full space M. d , the equivalence also 
allows us to define the hyperbolic Voronoi diagram on the boundary circle. Fig- 
ure 1 depicts the construction of the hyperbolic Voronoi diagram of a point set as 
the corresponding power diagram of associated balls. Note that the equivalence 
is only one way as there exists power diagrams with empty cells, which is not 
possible for the hyperbolic Voronoi diagrams. 



Hyperbolic Voronoi (Klein disk) Equivalent power diagram of balls 




Fig. 1. (Left) Hyperbolic Voronoi diagram in the Klein projective disk model is affine, 
and can thus be built equivalently from (right) a clipped power diagram of correspond- 
ing balls. Notice that the power diagram fully covers the Euclidean space. 

Building efficiently power diagrams is a well-studied problem [8] for which 
many highly optimized implementations are available. For example, CGAL 4 pro- 
vides this module. 

2.4 Converting Klein/Poincare disk models 

The hyperbolic Voronoi diagram in the Poincare disk P 2 is conformal, i.e. it pre- 
serves the notion of hyperbolic angles defined by intersecting curves in the rep- 
resentation space (the embedding space called the model). Straight lines in the 



4 The computational geometry algorithms library, see http://www.cgal.org 



Hyperbolic Voronoi diagrams made easy 5 



hyperbolic geometry are represented by arcs of circles perpendicular to the unit 
circle in the Poincare disk. Two non-intersecting arcs denote parallel hyperbolic 
lines, and two arcs intersecting orthogonally correspond to perpendicular lines. 
There exists infinitely many geodesies passing through a given point parallel to 
a given geodesic, contradicting Euclid's fifth postulate of Euclidean geometry. 
The distance between any two vector points p and q with Euclidean coordinates 
is h P (p,q) = arccosh(l + <5(p,q)) with <5(p,q) = 2 (1— | |p|f 2 T(i— | |q| - A geodesic 
(line segment in hyperbolic geometry) is expressed in the Poincare disk by either 
a portion of a circle x 2 + y 2 — 2(ax + by) + 1 = (with a 2 + b 2 > 1) or by a line 
ax = by. 

A point k in the Klein disk corresponds to a point p in the Poincare disk 
(and vice-versa) using the following radial scaling formula p = - — kf^ 'k anc ^ 



The relationship between the two Klein/Poincare models is therefore a pro- 
jection from the center of the disk. Figure 2 shows both Klein affinc and Poincare 
circular arc Voronoi diagrams. The straight line passing through two points in 
the Klein model intersects the unit circle in two ideal points, with the corre- 
sponding Poincare line being a circular arc orthogonal to the boundary disk. 



Fig. 2. Hyperbolic Voronoi diagram in the disk and upper plane models: (left) Poincare 
conformal Voronoi delimited by circular arcs superposed with the equivalent Klein 
affine Voronoi diagram, (b) Poincare upper plane conformal Voronoi diagram. 



2.5 Sphere inversions and Mobius transformations 

Mobius transformations (also called homographies) form a geometric group de- 
fined by inversions of spheres. An inversion of the space by a sphere centered at 
O with radius r, maps to itself all rays emanating from C such that the product 
of a point on that ray with its image (i.e., the mapped point) equals r 2 . Mobius 
transformations preserve angles as well. The isometries of hyperbolic geometries 



k = 
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are in one-to-one mapping with Mobius transformations. Thus we can navigate 
hyperbolic space by smoothly moving the viewpoint of the Poincare disk model 
using Mobius transformations (see Section 3.1). 

Note that in computational geometry, the Mobius Voronoi diagram [9] of 
a point set {p^}i is defined as the space partition relating proximity structure 
induced by distance functions -Dj(x) = Aj||x — pj 2 — fa, where the Aj's are 
multiplicative factor constants and the fa's are additive constants (i.e., additively 
and multiplicatively weighted Voronoi diagrams). Hyperbolic Voronoi diagrams 
are special cases of Mobius diagrams that can be viewed as affine diagrams. 

2.6 Converting to Poincare upper half plane 

To convert the affine Klein Voronoi diagram to the Poincare upper half plane 
model, we first convert the Klein disk model to the Poincare disk model as 
explained in the previous section. We then interpret 2D points (a, 6) as complex 
numbers z = a + ib in the complex space, so that the Poincare disk model is the 
space {zgC||2:|<l} and the upper half plane model corresponds to the space 
{z = a + ibEC\b> 0}. Both Poincare disk/upper plane models arc conformal. 
We then apply the following Mobius transformation / : C U {oo} — > C U {oo} 
with z i ► f(z) = ^jzr^ ■ Reciprocally, to convert the upper half plane to the disk 
model, we apply the inverse Mobius transformation / _1 : C U {oo} — > C U {oo} 

with z i ► r x {z) = §=j (with / o /- 1 = id). 

2.7 Hyperbolic Delaunay triangulations 

The power diagram equivalent to the hyperbolic Voronoi diagram admits a dual 
regular triangulation [8,2,10] as depicted in Figure 3. This straight line triangu- 
lation can be embedded in Klein disk model of projective geometry. Furthermore, 
the bisector of two points p and q is perpendicular to the geodesic linking p to q. 
In particular, the two circular arcs meet orthogonally in the Poincare conformal 
representations. In Klein disk model, they also meet perpendicularly but since 
the representation is not conformal, it makes a different Euclidean angle. How- 
ever, we can easily compute the angle between two line segments in Klein model 
by considering their chords. Indeed, the ideal points located on the boundary 
circle are the same in the Klein model and the Poincare disk model. 

2.8 Weighted and fe-order hyperbolic Voronoi diagrams 

Power diagrams are weighted (ordinary) Voronoi diagrams that can be obtained 
by projecting vertically the polytope defined by the intersection of half-spaces [8] 
generated by lifting points to the paraboloid z — (x,x) of R d+1 , taking the 
corresponding tangent plane and translating them vertically according to their 
weights. A similar lifting/projection mechanism holds for hyperbolic Voronoi 
diagram in the upper-plane model (the proof relying on pencils of spheres [11] 
is omitted) 
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Fig. 3. Hyperbolic Delaunay and regular triangulations obtained as dual structures of 
Voronoi diagrams. The regular triangulation (right) is embedded in the Klein model of 
hyperbolic geometry. 

Theorem 2. The 2D hyperbolic Voronoi diagram on the upper plane can be 
obtained by first computing the intersection of 3D halfspaces implied by the 
paraboloid lifting, then projecting horizontally the edges of the 3D polytope to 
the paraboloid, and finally projecting vertically these curves on the plane. 

Furthermore, as shown in [10], fc-order and weighted affinc diagrams are also 
affine diagrams that can therefore be computed equivalently as power diagrams. 
These power diagrams can be built themselves as the projection of a (d + 1)- 
dimcnsional polytope obtained as the intersection of half-spaces tangent to the 
paraboloid and translated vertically by a corresponding weight. 

3 Applications of the hyperbolic Voronoi diagram 

3.1 Browsing image collections in the hyperbolic disk 

Figure 4 shows a snapshot of our interactive image search and browser appli- 
cation. We use the fact the Poincare/Klein disk models allow one to display 
the full hyperbolic plane into a bounded disk fitting the display dimension. The 
application asks for text keywords and query web search engines to retrieve 
large collections of images. These images are analyzed using various feature ex- 
tractions (e.g., histograms) and represented in the 2D hyperbolic space using 
Sammon clustering technique 5 that optimizes the non-linear embedding so as to 
preserve pairwise distances [12] of high-dimensional extracted feature vectors. 

The Poincare/Klein disk models of the hyperbolic space are not uniquely 
defined as it depends on the chosen origin and orientation around the disk center: 
The selected viewpoint [13], also called focus [14]. From the user standpoint, it 
is important to be able to interactively change the focus so as to better visually 
depict the structures of the image collection. In particular, two useful user- 
interface operations involve computational geometry tasks 

5 See http://www.techfak.uni-bielefeld.de/~walter/h2vis/ 
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Fig. 4. Snapshot of the image browser applet in the Poincare hyperbolic disk. This 
bounded representation of hyperbolic space greatly helps in navigating among large 
image collections. Users enter keywords from which internet images are automatically 
aggregated via the Yahoo! API. The left and right snapshots are related by a Mobius 
transform driven by mouse user-interface interactions. 



— Image selection: Given a mouse click, report the closest anchored image, 

— Group selection: Given several mouse clicks, report the smallest enclosing 
ball containing all anchored images falling within the mouse selection. 

The first operation is a nearest neighbor query in the hyperbolic geometry. 
The second operation is useful for re-centering the hyperbolic disk model using 
the viewpoint so as to improve the display area of the selected images. Since 
rectangular shapes does not visualize rectangularly in either Poincare/Klein disk, 
we prefer to manipulate circles that are rasterized as circles in Poincare conformal 
disk and ellipsis in the Klein non-conformal disk. 

The transformations for re-centering the viewpoint of the disk models are 
Mobius transformations. 

Closest image and nearest neighbor queries. To answer the first specific 
UI operation of finding the nearest neighbor for a given mouse-click query, we 
build the Klein Voronoi diagram. This yields a 2D straight line partition of 
the disk space (i.e., an arrangement) on which we perform point location in 
logarithmic time [2]. We may recenter the viewpoint of the disk model at the 
selected image, if wished. 

Group selection and smallest enclosing ball. To address the second UI 
primitive, we first need to compute the smallest enclosing ball, and then recenter 
the viewpoint of the disk model so as to coincide with the circumcenter of that 
enclosing ball. The smallest enclosing ball may be computed from the furthest 
Voronoi diagram that is an affine diagram in Klein disk in 0(n log n) time, fol- 
lowing [15]. This is sub-optimal as we may rather use a simple linear randomized 
algorithm originally proposed by Welzl [16] for Euclidean geometry and recently 
extended to dually flat spaces [17]. Indeed the combinatorial optimization prob- 
lem is of type LP (linear programming). On the plane, we need to define the 
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two basic primitives: The hyperbolic ball passing through two points, and the 
hyperbolic ball passing through three points. The first case is solved by com- 
puting the intersection point of the Klein bisector with the line segment joining 
these two points. The second case considers the intersection any two of the three 
bisectors. Once the circumcenter of the smallest enclosing ball is computed, we 
finally recenter the viewpoint of the disk model using a corresponding Mobius 
transformation. 




Fig. 5. Smallest enclosing ball of a point set in the Klein non-conformal disk model: 
affine bisector and smallest enclosing ball deformed as an ellipsis. 



4 Conclusion 

We have shown that the hyperbolic Voronoi diagram in the Klein projective 
disk model is merely an affine diagram that can be conveniently computed as 
an equivalent power diagram. This crucial observation led us to a framework for 
performing other useful computational geometry tasks such as answering nearest 
neighbor queries or computing smallest enclosing balls. We illustrated how these 
techniques come in handy for designing user interfaces of an interactive image 
browser application. 
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